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Abstract: Using 328 pb^ of data collected at DA<I>NE corresponding to ~1.8 million 
Kl — > TTfii' decays, we have measured the form factor parameters. The structure 
of the K — TT vector-current provides information about the dynamics of the strong 
interaction; its knowledge is necessary for evaluation of the phase-space integral required 
for measuring the CKM matrix element Vus and for testing lepton universality in kaon 
decays. Using a new parametrization for the vector and scalar form factors, we find 
A+=(25.7 ± 0.6) X 10"^ and Ao=(14.0 ± 2.1) x lO'^. Our result for Aq, together with 
recent lattice calculations of /^r, Jk and /(O), satisfies the Callan-Treiman relation. 
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1. Introduction 




Semileptonic kaon decays, Kl 
(Fig. H) offer possibly the cleanest way to ob- 
tain an accurate value of the Cabibbo angle, or 
better, Vus- Since — > vr is a 0~ — >• 0~ transi- 
tion, only the vector part of the hadronic weak 
current has a non-vanishing contribution. Vec- 
tor transitions are protected by the Ademollo- 
Gatto theorem against SU(3) breaking correc- 
tions to lowest order in nis (or — mu,d)- At 
present, the largest uncertainty in calculating 
difficulties in computing the matrix element {■k\J]^^'^\K) . In the notation of Fig. |^ 
Lorentz invariance requires that this matrix element have the form 



Figure 1: Amplitude for Kl — > ir^i^v. 
The gray circle indicates the KttW ver- 
tex structure. 

from the decay rate is due to the 



TtI J^<^|i^) = iP+p)a f+{t) + (P - p)a f-{t) 



where t = [P — p)^ = {k + k')"^ = -|- m? — 2ME.,^ is the only L-invariant variable. 
The form factors (FF) f+{t) and f-{t) account for the non-pointlike structure of the 
hadrons; the values of the FFs at t = differ from unity because of SU{3) corrections, 
i.e., because pions and kaons have different structure. The term containing /_ is 
negligible for K^^ decays, because the coefficient P — p = k + k', when acting on the 
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leptonic current, gives the lepton mass. The FF /_ must be retained for i^^3 decays. 
It is customary to introduce a scalar FF /o(t), such that Eq. ( [L.l| ) becomes 



(^ijrii^> = /(o) 



with Ak-k = M'^ — m?. Since the FFs /+ and /o must have the same value at t = 0, 
the term /(O) has been factored out. The functions and /o(t) are therefore both 

unity at t = 0. If the FFs are expanded in powers of t up to as 

h,o{t) = 1 + A'+,o A + ^ ^+,0 (1-2) 

four parameters (A'_|_, A" , Aq, and Aq) need to be determined from the decay spectrum 
in order to be able to compute the phase space integral that appears in the formula 
for the partial decay width. However, this parametrization of the FFs is problematic, 
because the values for the As obtained from fits to the experimental decay spectrum 
are strongly correlated, as discussed in the Appendix and in Ref. [4]. In particular, the 
correlation between Ag and Ag is —99.96%; that between A^ and A^ is —97.6%. It is 
therefore impossible to obtain meaningful results using this parameterization. 
Form factors can also by described by a pole form: 

Mia 

f+{t),r,= , (1.3) 

which expands to 1 + t/My g + (t/Myg)'^, neglecting powers of t greater than 2. It is 
not clear however what vector and scalar states should be used. 

Recent K^^^ measurements [1-3] show that the vector FF is dominated by the closest 
vector (qq) state with one strange and one light quark (or Ktt resonance, in an older 
language). The pole- fit results are also consistent with predictions from a dispersive 
approach [5-7]. We will therefore use a parametrization for the vector FF based on a 
dispersion relation twice subtracted at t = [6]: 



f+{t) = exp 



-L{A+ + H{t)) 



:i.4) 



where H{t) is obtained using K — n scattering data, and /+(0) = 1, /:^(0) = A_|_/m^ . 
A good approximation to Eq. (|1.4D is 



~ , t Ai + 0.000584 / t y Al + 3 X 0.000584 A+ +0.0000299 / i 

U{t) = i+x„+^ — — +- + 



2 \m? J 6 \m? 

(1.5) 

The errors on the constants 0.000584 and 0.0000299 in Eq. (|]|) are 0.00009 and 
0.000002, respectively. The pion spectrum in decay has also been measured re- 
cently [1,8,9]. As discussed in the Appendix, there is no sensitivity to Ag. All authors 
have fitted their data using a linear scalar FF: 

/o(t) = l + Ao^. (1.6) 
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Because of the strong correlation between Aq and Aq, use of the hnear rather than the 
quadratic parameterization gives a value for Aq that is greater than Aq by an amount 
equal to about 3.5 times the value of Aq. To clarify this situation, it is necessary to 
obtain a form for fo{t) with t and terms but with only one parameter. 

The Callan-Treiman relation [10] fixes the value of scalar FF at t = Aktt (the so- 
called Callan-Treiman point) to the ratio of the pseudoscalar decay constants fx/fTr- 
This relation is slightly modified by SU(2)-breaking corrections [11]: 

fo{AKn) = yJ(q)+ ^ct, Act ^ -3.4 x lO'^ (1.7) 

A recent parametrization for the scalar FF [5] allows the constraint given by the Callan- 
Treiman relation to be exploited. It is a twice-subtracted representation of the FF at 
t = Aktt and t = 0: 

foit) = exp log(C - G(t))) (1.8) 

such that C = foiAxTr) and fo{t) = 1. G{t) is derived from Kn scattering data. As 
suggested in Ref. [5], a good approximation to Eq. (|1.8|) is 



t Ag + 0.000416 / t V A^ + 3 X 0.000416 Ao + 0.0000272 / t 
/Q(t) = l + Ao^ + (^^j + 

(1.9) 

with logC = XoAKTv/ml + 0.0398 ± 0.0041. Eq. (^) is quite similar to the result in 
Ref. [12]. The errors on the constants 0.000416 and 0.0000272 in Eq. (|L9|) are 0.00005 
and 0.000001, respectively. 

At KLOE, the pion energy and therefore t can be measured, since the momen- 
tum is known at a (j) factory. However, 7r-/i separation is very difficult at low energy. 
Attempts to distinguish pions and muons result in a loss of events of more than 50% 
and introduce severe systematic uncertainties. We therefore use the neutrino spectrum, 
which can be obtained without vr-// identification. 



2. The KLOE detector 

The KLOE detector consists of a large cylindrical drift chamber, surrounded by a 
lead scintillating-fiber electromagnetic calorimeter. A superconducting coil around the 
calorimeter provides a 0.52 T field. The drift chamber [13] is 4 m in diameter and 
3.3 m long. The momentum resolution is (Tp^/p± ~ 0.4%. Two-track vertices are 
reconstructed with a spatial resolution of -^3 mm. 

The calorimeter [14] is divided into a barrel and two endcaps. It covers 98% of 
the solid angle. Hits on cells nearby in time and space are grouped into calorimeter 
clusters. The energy and time resolutions are oe/E = 5.7%/y^ E (GeV) and ctt = 
54 ps/y^S (GeV) © 100 ps, respectively. 

The KLOE trigger [15] uses calorimeter and chamber information. For this analysis, 
only calorimeter information is used. Two energy deposits above threshold {E > 50 
MeV for the barrel and E > 150 MeV for endcaps) are required. Recognition and 
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rejection of cosmic-ray events is also performed at the trigger level. Events with two 
energy deposits above a 30 MeV threshold in the outermost calorimeter plane are 
rejected. 

3. Analysis 

Candidate Kl events are tagged by the presence of a Ks — > tt'^tt^ decay. The 
tagging algorithm is fully described in Refs. 16. The momentum, pkl^ is obtained 
from the kinematics of the (j) KsK^ decay, using the reconstructed Ks direction and 
the known value of p^. The resolution is dominated by the beam-energy spread, and 
amounts to about 0.8 MeV/c. The position of the (j) production point, x^, is taken as 
the point of closest approach of the Kg path to the beam line. The Kl line of flight 
{tagging line) is given by the Kl momentum, p;^^ — PXg and the position of the 
production point, x,^. All tracks in the chamber, after removal of those from the Kg 
decay and their descendants, are extrapolated to their points of closest approach to the 
tagging line. 

For each track candidate, we evaluate the point of closest approach to the tagging 
line, Xc, and the distance of closest approach, dc- The momentum pc of the track at Xc 
and the extrapolation length, Ic, are also computed. Tracks satisfying dc < ar^y + b, 
with a = 0.03 and 6 = 3 cm, and —20 < /c < 25 cm are accepted as Kl decay products. 
Vxy is the distance of the vertex from the beam line. For each sign of charge, we chose 
the track with the smallest value of dc as a decay product, and from them we 
reconstruct the decay vertex. Events are retained if the vertex is in the fiducial volume 
35 < r^y < 150 cm and \z\ < 120 cm. The combined tracking and vertexing efficiency 
for Kfj_2 is about 54%. This value is determined from Monte Carlo (MC), corrected 
with the ratio of data and MC efficiencies obtained from Kl—>-tt^it~tt^ ,7rei^ control 
samples [16]. 

Background from Kl^tt^'k~ , 7r+7r~7r'' is easily removed by loose kinematic cuts. 
The largest background is due to Kl ir^e^u decays, possibly followed by early vr — > 
HI' decay in flight. For all candidate i^^s events we compute min(A7re, Aett), the smaller 
value of l-Emiss — PmissI assumiug the decay particles are vre or en. We retain events only 
if this variable is greater than 10 MeV. After the above kinematic cuts the efficiency 
for the signal is about 96% and the purity is about 80%. 

A further cut on the scatter plot of A^t^ = E^issifJ-'^ ,'^~) — Pmiss vs A^^^ = 
-Emiss(''r'^, — Pmiss shown in the left panel of Fig. ^ for Kl nfiu and background 
events respectively, is applied. The right panel of Fig. ^ shows the distribution of 
the lesser between A,r/i and A^^^ for data and MC. After the kinematic cuts described 
above, the contamination, dominated by Ki^nev decays is ~4%. 

To further reduce Ki^nev background we use the particle identification (PID) 
based on calorimeter information. Tracks are required to be associated with EMC 
clusters. We define two variables: dTC) the distance from the extrapolated track entry 
point in the calorimeter to the cluster centroid and d^^TC) the component of this 
distance in the plane orthogonal to the track momentum at the calorimeter entry point. 
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Figure 2: Left: A^^r versus A^r^ distribution from MC. Kj^ -^t^\xv (gray scale) and background 
(black points). The outsermost contour shows the accepted region. Right: r7im(A^^, A^^) for 
data (black dots), MC (solid line), and MC signal (gray shaded histogram). 



: DATA 
e MC 



We accept tracks with d^^TC < 30 cm. The cluster efficiency is obtained from 
the MC, corrected with the ratio of data and MC efficiencies obtained from control 
samples. These samples, of 86% and 99.5% purity, are obtained from ET^s and iCes 
events selected by means of kinematics and independent calorimeter information. The 
cluster efficiency correction versus Ey is shown in Fig. |^ 

For each Ki^ decay track with an associated 
cluster we define the variable: Atj = tci — ^i-, {i = 
TT, e) in which tc\ is the cluster time and ti is the 
expected time of flight, evaluated using the corre- 
sponding mass, ti includes the time from the entry 
point to the cluster centroid [17]. We determine 
the e^e~ collision time, to, using the clusters from 
the Ks. 

The mass assignment, vre or en, is obtained by 
choosing the lesser of |At^+ — Atg-| and lAt^^- — 
Atg+ 1. After the mass assignment has been made, 
we consider the variable 
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Figure 3: Cluster efficiency correc- 
tion versus E^. 

At^ - Atr 
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where 0"+ (o"_) = 0.5(0.4) ns are the resolutions. Additional information is provided by 
the energy deposition in the calorimeter and the cluster centroid depth. These quanti- 
ties are input to a neural network (NN). We retain events with VFmax < (1/6) Rtof+0-^, 
where Wmax is the largest of the NN outputs for the two charge hypothesis. The dis- 
tribution of Wmax -l/6i?TOF — 0.4 for data and MC is shown in Fig. |^. The resulting 
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Figure 4: Distribution of Wmax -l/6i?TOF — 0.4 for data (dots), MC (solid line) and MC signal 
(gray scale). The dashed line indicates the cut that we use. 



purity of the sample is ~ 97.5%, almost uniform in range 16< E^y <181 MeV used for 
the fit. 

The FF parameters are obtained by fitting the E^^ distribution of the selected 
events in the range 16 < E^ < 181 MeV, subdivided in 32 equal width bins. The 
bin size, 5 MeV, is about 1.7 times the neutrino energy resolution. The value of E,^, 
i.e. the missing momentum in the Kl rest frame, is determined with a resolution of 
about 3 MeV almost independently on its value. The purity of the final sample used 
to extract the form factor parameters is illustrated in Fig. |5|. 

After subtracting the residual background 
as estimated from MC, we perform a fit to 
the data using the following expression for the 
expected number of events in each of the 32 bins: 



1.00 



(3.1) 




20 40 60 80 100 120 140 160 180 200 
Figure 5: Purity versus E^. 



where Arj(A) is the fraction of events expected 
for the parameter set defined by A in the j*^ bin, 
and Aij is the resolution smearing matrix. FpsR 

is the final state radiation correction. It is evaluated using the MC simulation, GEANFI 
[18], where radiative processes are simulated according the procedures described in 
Ref. [19]. FSR affects E'l.-distribution mainly for high E^y values, where the correction 
is about 2%. The free parameters in the fit are the FFs A. No , the total number of 
signal events, is fixed. 
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4. Systematic uncertainties 

The systematic errors due to the evaluation of corrections, data-MC inconsistencies, 
result stability, momentum mis-calibration, and background contamination are sum- 
marized in Tab. [l|, for the case of a quadratic and a linear fo{t) ■ 



Source 


X 10^ 


6X1 X 


SXo X lO-'^ 


Tracking 


1.60 


0.47 


0.86 


Clustering 


2.07 


0.61 


1.87 


TOF + NN 


2.23 


1.16 


1.45 


p-scale 


1.10 


0.71 


0.81 


p-resolution 


0.61 


0.21 


0.01 


Total 


3.66 


1.58 


2.64 



Table 1: Summary of systematic uncertainties on A'+, A" , Aq. 



The uncertainty on the tracking efficiency correction is dominated by sample statis- 
tics and by the variation of the results observed using different criteria to identify tracks 
from Kl decays. Its statistical error is taken into account in the fit. We study the effect 
of differences in the resolution with which the variable dc is reconstructed in data and in 
MC, and the possible bias introduced in the selection of the control sample, by varying 
the values of the cuts made on this variable when associating tracks to Kl vertexes. For 
each variation, corresponding to a maximal change of the tracking efficiency of about 
±10%, we evaluate the complete tracking-efficiency correction and measure the slope 
parameters. We observe changes of 1.60x10^^, 0.47x10"^, and 0.86x10"^ for A'^ A+ 
and Ao , respectively. 

As for tracking, we evaluate the systematic uncertainties on the clustering efficiency 
corrections by checking stability of the result when the track-to-cluster association 
criteria are modified. The statistical uncertainty on the clustering efficiency corrections 
is taken into account in the fit. The most effective variable in the definition of track- 
to-cluster association is the transverse distance, d^^TC- We vary the cut on (i_L^TC 
in a wide range from 15 cm to 100 cm, corresponding to a change in efficiency of 
about 19%. For each value of the cut, we obtain the complete track extrapolation and 
clustering efficiency correction and we use it to evaluate the slopes. We observe changes 
of 2.07x10-3, 0.61x10-3 and 1.87x10-3 for A'+ A'| and Aq, respectively. 

We study the uncertainties of the efficiency and of the background evaluation by 
studying the stability of the result with modified PID and kinematic cut values, cor- 
responding to a variation of the cut efficiency from 90% to 95%. This changes the 
background contamination from 1.5% to 4.5%. We observe changes of 2.23x10^^^ 
1.16x10-3 and 1.45x10-3 for A'^_ A'^ and Aq, respectively. 

We also consider the effects of uncertainties in the absolute momentum scale and 
momentum resolution. A momentum scale uncertainty of 0.1% [18] corresponds to 
changes of 1.1 x 10-3, 0.71 x 10-3, and 0.81 x 10-3 for A'_,_, A'|, and Aq, respectively. 
We investigate momentum resolution effects by changing the value of the resolution 
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on E,y by 0.1 MeV, an amount which characterizes our knowledge of the momentum 
resolution, as described in Ref. [3]. We observe changes of 0.61 x 10~^, 0.21 x 10~^, 
and 0.01 x 10^'^ for A', , A'l , and Aq, respectively. 



5. Results and interpretation 



About 1.8 million decays were accepted. We first fit the data using Eqs. (1.2) 
and ( |1.6D for the vector and scalar FFs. The result of this fit is shown in Fig. ^. 
We obtain: 

0.02 


-0.02 



AV = (22.3±9.8stat±3.7syst) x 10 



A^; = (4.8±4.9stat±1.6syst) 
Ao = (9.1±5.9stat±2.6syst) 



10^ 
10" 



-3 
3 
3 




80000 - 
60000 - 
40000 - 
20000 - 



(Data-Fit) /Fit 



Events 



Fit 
Data 



(MeV) 



with x^/dof = 19/29, and correlation 
coefficients as given in the matrix. 

Improved accuracy is obtained by 
combining the above results with those 
from our K^^ analysis [3]: 

AV = (25.5 ± 1.5stat ± l.Osyst) X 10-3 
A^; = (1.4±0.7stat ±0.4syst) X 10-3 



20 40 60 80 100 120 140 160 180 200 

Figure 6: Residuals of the fit (top plot) and Ei, 
distribution for data events superimposed on the 
fit result (bottom plot). 



We then find: 

AV = (25.6 ± 1.5stat ± 0.9syst) X lO'^ 
A'l = (1.5 ± 0.7stat ± 0.4syst) X 10-3 
Ao = (15.4 ± l.Sstat ± 1.3syst) X 10-3 

with x^/dof = 2.3/2 and the correlations given in the matrix on the right. 

Finally, to take advantage of the recent parameterizations of the FFs based on 
dispersive representations (Eqs. (|1.4D and ( |1.8[ )), we combine our results from this 
analysis of K^^ data with our previous result for Kgs. We perform a fit to the values 
obtained for A'_,_, A'^, and Aq that makes use of the total error matrix as described above, 
and the constraints implied by Eqs. ( |l.5D and (|l.9[) . Thus, the vector and scalar FFs 
are each described by a single parameter. Dropping the " ' " notations, we find 

A+ = (25.7 ± 0.4stat ± 0.4syst ± 0.2param) X 10-3 
Ao = (14.0 ± 1.6stat ± 1.3syst ± 0.2param) X 10-3 
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with x^/dof = 2.6/3 and a total correlation cofficient of —0.26. The uncertainties 
arising from the choice of parameterization for the vector and scalar FFs are 0.2 x 10~^ 
and 0.1 X 10^'^ using only ivTgs decays and K^^ decays, respectively. These contributions 
to the uncertainty on and Aq are explictly given in Eq. dSTl). We note that the use 
of Eq. ( |1.9| ) changes the value of the phase space integral by only 0.04% with respect 
to the result obtained using a linear parameterization for /o(t)- 

Finally, from the Callan-Treiman relation we compute /(O) = 0.964 it 0.023 using 
fx/ fn = 1-189 lb 0.007 from Ref. [20]. Our value for /(O) is in agreement with the 
results of recent lattice calculations [21]. 



6. Conclusions 

We have performed a new measurement of the Kl — > vr/xi^ FFs. Our results are in 
acceptable agreement with the measurements from KTeV [1] and ISTRA+ [9] and in 
disagreement with those from NA48 [2,8]. In particular, our result for the scalar FF 
parameter Aq = 0.0143 it 0.00203 (Eq. accounts for the presence of a term. 

KTeV and ISTRA+ use a linear parameterization; as a consequence, their values for 
Ao are systematically high by ~0.003. We also derive f^\o) = 0.964 ± 0.023. This 
value is in agreement with the results of recent lattice calculations [21]. 



A. Error estimates 



It is quite easy to estimate the ideal error in the measurements of a set of parameters 
p=(pi, p2, ■ ■ ■ Pn) from fitting some distribution function to experimentally determined 
spectra. Let F{p, x) be a probability density function, PDF, where p is some parameter 
vector, which we want to determine and x is a running variable, like t. The inverse of 
the covariance matrix for the maximum likelihood estimate of the parameters is given 
by [22]: 

(G-i),. = -^!M 
dpidpj 

from which, for N events, it trivially follows: 

(G~^).. = N — — — dv, 
V J F dpi dpj ' 

with dv the appropriate volume element. We use in the following the above relation 
to estimate the errors on the FF parameters for one and two parameters expression 
of the FFs and /o(t)- The errors in any realistic experiment will be larger than 

our estimates, typically two to three times. The above estimates are useful for the 
understanding of the problems in the determination of the parameters in question. 
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A.l Ke3 decays 

For a quadratic FF, f{t) = l + X'_^{t/m'^) + {X'^/2){t/m'^)^, the inverse of the covariance 
matrix GT^, the covariance matrix G4. and the correlation matrix are: 



N 



5.937 13.867 \ 
13.867 36.2405 / 



1.258^ -0.606" 
-0.606 0.5092 



1 




The square root of the diagonal elements of G+ gives the errors, which for one milhon 
events are (5AY=0. 00126, =0.00051. The correlation is very close to —1, meaning 
that, because of statistical fluctuation of the bin counts, a fit will trade for A" and 
that the errors are enlarged. A fit for a linear FF, f(t) = 1 + X'_^{t/m?) in fact gives 
A^=0.029 instead of 0.025 and an error smaller by ~3: 



5AV = VG+(1,1) = 0.0004. 

A simple rule of thumb is that ignoring a t'^ term, increases A^ by ~3.5xA" . For K^s 
decays the presence of a t"^ term in the FF is firmly established. It is however not fully 
justified to fit for two parameters connected by the simple relation X'^=2xX'_^'^ . The 
authors of ref. 5, 6 explicitly give an error for their estimate of the coefficient of the 



t terms. The above discussion justifies the use of eq. 1.5. The errors obtained above 



compare reasonably with the errors quoted in [1-3], when all experimental problems 
are taken into account. 

A. 2 K^3 decays 

The scalar FF only contributes to K^^ decays. Dealing with these decays is much harder 
because: a) - the branching ratio is smaller, resulting in reduced statistics, b) - the E.,^ 
or t range in the decay is smaller, c) - it is in general harder to obtain an undistorted 
spectrum and d) - more parameters are necessary. This is quite well evidenced by the 
wide range of answers obtained by different experiments [1,8,9]. Assuming that both 
scalar and vector FF are given by quadratic polynomials as in eq. ordering the 
parameters as Ag, Ag, A'^_ and A" , the matrices Gq^_j^ and Go& + , are: 



N 



and the correlations, ignoring the diagonal terms, are: 

/ -0.9996 -0.974 0.91 \ 
0.978 -0.919 
-0.976/ 



/1.64 


5.44 


1.01 


3.90 \ 




/ 63.92 


-1200 


-923 


197 \ 


5.44 


18.2 


3.01 


12.3 


1 


-1200 


18.82 


272 


-59 


1.01 


3.01 


1.47 


4.24 


' N 


-923 


272 


14.82 


-49 


\3.90 


12.3 


4.24 


13.8 / 




V 197 


-59 


-48 


3.42/ 



V 



(A.l) 



All correlations are very close to —1. In particular the correlations between Ag and Aq 
is —99.96%, reflecting in vary large 5X'q and 5Xq errors. We might ask what the error on 
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Aq and Aq might be if we had perfect knowledge of A^ and A^. The inverse covariance 
matrix is give by the elements (1,1), (1,2), (2,1) and (2,2) of the Gq^_(_ matrix above. 
The covariance matrix therefore is : 



Go(A^, A^' for A'+, A'| known) 



For one million events we have (5Aq=0.0024, about 4x the expected value of Aq. In 
other words Aq is likely to be never measurable. It is however a mistake to assume 
a scalar FF linear in t, because the coefficient of t will absorb the coefficient of a 
term, again multiplied by '^3.5. Thus a real value Aq=0.014 is shifted by the fit to 
0.017, having used eq. 1.9. Fitting the pion spectrum from 1 million Kf^^ decays for 
Ao, A+ with the FFs of eq. |l| and |l| gives the errors 5Ao~0.0096 and (5AY~0. 00097. 
Combining with the result from a fit to 1 million K^^, with the FF of eq. L5 for 
which (5A'+'-0.00037 gives finally (5A'o~0.00075, (5A'+~0. 00034 and a Ao-A'+ correlation of 
—31%. Using the neutrino spectrum for Kf^^ decays, the errors are only slightly larger: 



(5Aq~0.001, 5A'_)_~0. 00036. We hope to reach this accuracy with our entire data sample, 
~5x the present one, and a better analysis which would allow using the pion spectrum. 
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